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Stochastic amplification of weak signals in an RF SQUID with ScS contact 
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The paper presents results of a numerical simulation of the stochastic dynamics of magnetic flux in an RF 
SQUID loop with a Josephson point (ScS) contact driven by a mix of band-limited Gaussian noise and 
low-frequency small-amplitude sine signal, at finite temperatures 0 < T < T c . A change in the gain of the 
stochastically amplified weak signal with the temperature rise is shown to be due to smearing of the energy 
barrier between adjacent metastable current states of the loop. A comparison to an RF SQUID with a tunnel 
(SIS) junction is done. 
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I. INTRODUCTION 

The magnetometers based on Superconducting Quan¬ 
tum Interference Devices (SQUIDs) are widely used in 
physical experiments, medicine (magnetocardiographs 
and magnetoencephalographs), geophysics (underground 
radars), electronic manufacturing (SQUID microscopes). 
The sensitivity of SQUIDs and their quantum ana¬ 
logues, SQUBIDs, has practically reached the quantum 
limitation^ -1. However, with increase of the quantizing 
loop inductance up to L ~ 10 -9 — 10“ 10 H , the ther¬ 
modynamic fluctuations lead to expressed degradation of 
the energy resolution. 

As shown earlier—the sensitivity of the magnetome¬ 
ters can be considerably enhanced in this case by us¬ 
ing the stochastic resonance (SR) effect. The SR phe¬ 
nomenon discovered in the early 1980s^ manifests itself 
in a non-monotonic rise of a nonlinear (often bi-stable) 
system response to a weak periodic signal which peaks 
at a certain intensity of the noise added (or inherent) to 
the system. 

Owing to extensive studies during the last two decades, 
the stochastic resonance effect has been revealed in a 
variety of natural and artificial systems, both classical 
and quantum ones. Analytical approaches and quan¬ 
tifying criteria for estimation of the noise-induced or¬ 
dering were determined and then summarized in the 
reviewaiSiii. For example, the gain of 40 dB was ex¬ 
perimentally demonstrated^ for a weak harmonic signal 
due to SR in a SQUID with the tunnel (SIS) Josephson 
junction. For an under-optimal noise intensity, the sig¬ 
nal gain in SR SQUIDs can be nevertheless maximized 
with the stochastic-parametric resonance (SPR) effect— 
emerging when the system is affected by the weak signal, 
the noise and an additional high-frequency electromag¬ 
netic field simultaneously. 

Not so long ago, researchers paid their attention 
to clean small-sized superconducting contacts of ”con- 
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striction” (ScS) type called the ’’atomic-size contacts”, 
ASCsii They have only a few quantum conductive chan¬ 
nels in their cross section so they often referred to as the 
’’quantum point contacts”, QPCs 14 . The critical current 
I c of such contacts may take discrete values. The interest 
to the QPCs is roused by studies of the channel quan¬ 
tum conductivity and building superconducting qubits 
with large energy level splitting A E/h ~ 30 GHz£. The 
current-phase relation (the supercurrent I s as a function 
of the order parameter phase difference ip ) in the clean 
ScS contacts with ballistic electron transit at low tem¬ 
peratures ( T —> 0 ) 14 : 15 substantially differs from the 
classical Josephson current-phase relation for a SIS junc¬ 
tion I^ IS (p) = I c siny: and features a ’’saw-like” shape 
with discontinuities at p = mr. 

The coupling energies for SIS and ScS obtained by in¬ 
tegrating the current-phase relations therefore differ, too. 
If the SIS junction is incorporated into a superconducting 
loop and an external magnetic flux 4> e = 4>o/2 is applied 
( <f>o = h/2e S3 2.07• 10 -15 Wb is the magnetic flux quan¬ 
tum), its current-phase relation forms the symmetric two- 
well potential of the loop, which is needed for SR, only 
at sufficiently large loop inductances and/or the junc¬ 
tion critical current, namely, when (3l = 27rL/ c /<f>o > 1 . 
Unlikely, the potential energy U 3cS (<&) of the supercon¬ 
ducting loop with a QPC has a ’’sharp” barrier with the 
singularity in its top separated two metastable current- 
carrying states of the loop with different intrinsic mag¬ 
netic flux $ for any /3l including /3l « 1- 

In classical case for the zero-temperature limit, the 
stochastic dynamics of the magnetic flux in a super¬ 
conducting loop with the Josephson ScS contact (ScS 
SQUID) substantially differs from that for SIS SQUID 
studied earlier^i Taking into account the quan¬ 
tum fluctuations at finite temperatures 0 < T < T c 
changes the current-phase relation for the ScS contact 
and smoothes the potential barrier in the loop with ScS 
contact—. In this work we analyze the stochastic amplifi¬ 
cation of weak low-frequency harmonic signals in a super¬ 
conducting loop interrupted by Josephson ScS contact, 
at finite temperatures and various /3l- The results are 
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compared with the stochastic dynamics of a SIS SQUID. 


II. MODEL 

The stochastic dynamics of the magnetic flux inside an 
RF SQUID loop (inset in Fig.[T|) is studied by numerical 
solving the motion equation (Langevin equation) in the 
model of resistively shunted junction (RSJ model): 



x=<t >/® 0 

FIG. 1. The potential energy of RF SQUID loop with ScS 
contact (solid line) and SIS junction (dotted line) as a func¬ 
tion of reduced magnetic flux in the loop. The non-linearity 
parameter /3l = 12, fixed external magnetic flux = $o/2. 
The inset shows schematically the RF SQUID loop with a 
Josephson junction, the designations are described in the text. 
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has the ” saw-like” shape with discontinuities at ip = mr. 
Correspondingly, the potential energy of the supercon¬ 
ducting loop with ScS contact is 
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where Ej cS = J c $ 0 / 7 r is the maximal coupling energy of 
the Josephson ScS contact. U ScS ($, <f> e ) has singulari¬ 
ties in tops of the barriers separating the adjacent local 
minima which correspond to the loop metastable current 
states. 

Introducing the dimensionless parameter of non¬ 
linearity 


Pl = 2nLI c /<S> 0 , (5) 

and reducing the fluxes by the flux quantum $o : % = 
<!>/<I>o , x e = $e/4 ) o, and the potential energy by $g/2 L, 
expressions m and a can be rewritten, correspondingly, 
as 
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where C is the capacitance, Rn and I c are the normal 
shunt resistance and the critical current of the Joseph¬ 
son junction, respectively, $(t) is the magnetic flux inside 
loop, U (d>, $ e ) is the loop potential energy of the which is 
the sum of the loop magnetic energy Um and the Joseph¬ 
son junction coupling energy Uj. The time-dependent 
external magnetic flux <1 > e (t) applied to the loop contains 
the fixed and the varying components including the noise 
one. The junction coupling energy Uj is type-determined 
(for SIS, ScS, SNS, etc.). The potential energy of a loop 
with the ’’traditional” tunnel (SIS) junction is equal toi£ 
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where Ej IS = I c $o/2tt is maximal coupling energy of 
the tunnel Josephson junction. 

We consider the dynamics of the magnetic flux in the 
RF SQUID loop incorporating the clean ScS contact with 
ballistic electron transit*^. 

For classical^ and quantum^ ScS contacts with the 
critical current I c the current-phase relation in the zero- 
temperature limit T = 0 


SrS / \ ^e) pL I I /rr\ 

U ScS { X,X e ) = ---— IcOSTTaj. (7) 

The potential energy u SIS (x,x e ) of the loop with a 
tunnel junction has two or more local minima only at 
Pl > 1- It can be symmetrized by applying fixed mag¬ 
netic field $ e = ( x e = 1/2) . Fig.[T]shows potential 

energy of the loops with SIS and ScS contacts as a func¬ 
tion of internal magnetic field at large Pl = 12 for better 
visibility. The distinctive feature of the potential energy 
u ScS (x, x e ) of an RF SQUID loop with the ScS contact is 
the finite height of the ’’sharp” barrier between adjacent 
states for any vanishingly small Pl, he. L and/or I c (at 
zero temperature). Fig. [2] displays potential energies of 
the RF SQUIDs with SIS or ScS contact at several small 
Pl- It is obvious from Fig. [2] that the barrier separat¬ 
ing the two metastable current-carrying states vanishes 
at Pl = 1 in the SIS SQUID while stays finite in the ScS 
SQUID. 

The noise of thermal or any other origin causes switch¬ 
ing between the metastable states corresponding to the 
minima of U( ( h) . The switching rate for the white Gaus¬ 
sian noise with intensity (variance) D for parabolic wells 
and A U/D » 1 is given by the Kramers’ formula^: 


VK = tq exp(—AI//.D). 


( 8 ) 
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FIG. 2. The potential energy of the ScS (solid line) and SIS 
(dashed line) SQUIDs as a function of internal reduced mag¬ 
netic flux for several /3l ■ The fixed external flux 4> e = 4*o/2 . 
The barrier height AU and the inter-minima spacing Aa: are 
shown in the curves for (3l = 1.5 . 

In the case of thermal noise for T = 0 , D = 2fcgT. 
In this work we do not focus on any specific nature of 
the noise. The only requirement imposed to ensure the 
adiabatic mode of the SQUID switching is the limitation 
of the noise frequency band by a frequency not exceeding 
the reciprocal relaxation time for the internal flux in the 
SQUID loop 1/tl = R/L. Although the pre-exponential 
factor 7'o is different for the ” smooth” potential of the SIS 
SQUID and ’’sharp” potential of the ScS SQUID (details 
can be found in^), the exponent factor dominates any¬ 
way. Adding a small periodic signal with frequency f s to 
the external flux <f* e results in the stochastic-resonance 
system motion in the noise-driven bistable potential at 


and with 0 for the ScS contact, as 

dx 1 ... 0 L . . . . . . 

— = — 1 x e (t) — x + —(bvo.'kx) sgn(cos'KX)\. (12) 

dt T "[j 27T 

The external flux x e (t) is the sum of the fixed bias 
flux Xdc = 0.5, the useful signal x ac = asin2nf s t and 
the noise flux Xn- In the theory, the noise is often 
treated as the cj-correlated white, Gaussian-distributed 
noise: Xn = £(f), (£(t) ■ £(i — t')) = 2DS(t — t'). Dur¬ 
ing the numerical simulation, it is modeled by the quasi¬ 
random number generator giving Gaussian-distributed 
values with variance D = a 1 and the repetition period 
of about 2 90 . The sampling rate is 2 16 when solving the 
equation by finite differences technique that corresponds 
to the equivalent noise frequency band 32 kHz. So, the 
noise can be considered as quasi-white for the stochastic 
amplification of the signal with frequency f s = 10 Hz. 

Eqs. (fill) and (fl2l) were solved by the Heun method 
modified for the stochastic equations^. Each run gave 
an 8-s-long time series with its individual noise realiza¬ 
tion and then were subjected to fast Fourier transform 
(FFT). The obtained output spectral densities S<s>( oj) of 
the internal flux were averaged over 30 runs. The spec¬ 
tral amplitude gain of the weak periodic signal was de¬ 
termined as the ration of spectral densities of the input 
(external) and output (internal) magnetic fluxes: 

=syiH/syiM. (is) 

III. RESULTS AND DISCUSSION 


r K ~ 2 f s . ( 9 ) 

Taking typical experimental figures L = 3 • 10 _1 ° H, 
C = 3 • 10" 15 F, R n = 1...10 2 n, I c = 10 -5 ...10 -6 A, 

we found the McCamber parameter describing the effect 
of capacitance /3c = 2 ttR 2 I c C/$o < 1. This case corre¬ 
sponds to the viscous, non-oscillatory, motion and there¬ 
fore the term with the second derivative in Eq. © can 
be neglected. Low signal frequency f s = 10 Hz « 1/tl 
and noise cut-off frequency f c ~ 10 4 Hz « 1 /tl makes 
the problem adiabatic. This allows us the time depen¬ 
dence of the external flux consider as the time-dependent 
potential energy thus writing the motion equation in the 
form 


The current-phase relation for the ScS contact at a 
non-zero temperature T 14 ’ 15 


/f cS ((p) = I c S in | tanh 


A (T)cos f 
2 k B T 


(14) 


where lf cS {ip) is the contact supercurrent, I c = is 

the contact critical current, A(T) is the superconductor 
energy gap (order parameter), is the order parameter 
phase difference over the contact, ks is the Boltzmann 
constant, e is the electron charge, Rn is the contact nor¬ 
mal resistance, smoothes, tending, as approaching the 
superconductor critical temperature T c , to the sine law 
that describes the SIS junction. Consequently, the (re¬ 
duced) potential energy of the loop with ScS contact 


dx du(x , t) 

TL Tt + —&r = Q - 
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With (0 for the SIS junction, the Eq. m reads as 


= — [x e (t) - x + ^-sm2nx], (11) 

dt tl 47t 


„ ?S = (fLXfA ( 2 c« h ^), (15) 

where t = - T^ Tc , also becomes closer in its 

1.76-/1 — (T/T c ) 4 

shape to the SIS SQUID potential energy. The singu¬ 
larity in the barrier top disappears, while the barrier 
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height decreases with the temperature rise. Fig. [3] illus¬ 
trates the temperature evolution of ScS SQUID potential 
energy. The SIS SQUID potential energy is also given 
in the Figure for the comparison sake. In contrary to 
the zero-temperature case, the potential barrier vanishes 
at a certain reduced temperature even for (3l > 1 (at 
T/T c = 0.5 for the chosen /3 = 1.2) The inter-minima 
spacing decreases, too. The barrier heights for the ScS- 
and SIS-SQUIDs become equal at a certain temperature 
(at T/T c = 0.375 for the given case). 



FIG. 3. The potential energy of ScS SQUID as a function of 
the magnetic flux in the loop at various reduced temperatures. 
Similar function for the SIS SQUID at the same /3l is shown 
for comparison. 

By substituting (fl5l) in m, we get the motion equa¬ 
tion for the reduced flux in ScS SQUID at finite temper¬ 
atures T , 0 < T < T c , similar to Eq. m for the case 
T = 0. The solution of this equation at various noise 
intensities with further FFT as described above gives the 
spectral gain of a weak signal as a function of the noise 
intensity, i.e. the SR curves. 

Fig. [4] shows the SR curves in the ScS SQUID calcu¬ 
lated for several reduced temperatures and /3 l = 1-21. 
Strictly saying, /3l parameter is indicated for T = 0, since 
it renormalizes as the temperature rises. For the compar¬ 
ison, the SR curve in a SIS SQUID loop is displayed in the 
figure for the same /3l- One can see that, with the tem¬ 
perature rise, the maximal gain in ScS SQUID increases 
due to lowering of the barrier. This is much like a SIS 
SQUID behavior when /3l reduces down to the unityi&, 
so we consider the situation as a ’’renormalization” of /3l 
in ScS SQUIDs. 

When the barriers in the ScS and SIS SQUIDs are of 
the same height ( T/T c = 0.375, see Fig.[3|), the maxi¬ 
mal gain is observed at equal noise intensities for both 
SQUIDs while the maximal gain values are different be¬ 
cause of difference in the potential shapes and inter¬ 
minima spacings 

According to the two-state theory for linear response^ 
which does not accounts for inter-well dynamics, the 
weak-signal gain is determined by the spacing Ax be¬ 
tween the adjacent local minima of the system potential 
energy. Taking the SR curves in Fig.Q] for ScS and SIS 



Noise intensity s 


FIG. 4. The weak harmonic signal gain as a function the 
noise intensity in SIS and ScS SQUIDs at various reduced 
temperatures. The signal amplitude is 0.001, frequency 10 

Hz. 


SQUIDs with almost equal maximal gain ( T/T c = 0.45 
for ScS SQUID) and looking up corresponding potential 
energies in Fig.[4l we can see that Ax values are some¬ 
what different. This implies that the exact shapes of 
the potential profile should be involved in accurate the¬ 
oretical calculation. Nevertheless, it can be generally 
stated that the stochastic amplification of an informa¬ 
tional signal in ScS SQUIDs with small will reveal 
its specific features at ultra-low temperatures. The cal¬ 
culation shows the potential barrier in ScS SQUID with 
/3l = 0.1 vanishes at T/T c = 0.045 , that for niobium 
with its critical temperature T c = 9.25 K corresponds to 
a millikclvin temperature range. At such low tempera¬ 
tures, quantum effects will give a dominant contribution 
in switching rate22. Moreover, the ’’sharp” potential of a 
superconducting loop with ScS contact makes quantum 
tunneling between separated wells much faster as com¬ 
pared to a ’’common” SIS circuit. This feature allowed 
authors of£ to propose a superconducting qutrit-detector 
with tunneling rate of 20-30 GHz (which was important 
there for large energy level splitting). Thus, an interest¬ 
ing situation could emerge in ScS SQUIDs at ultra-low 
temperatures when two characteristic time scales would 
occur in the context of stochastic resonance, one for clas¬ 
sical hops over the barrier and other for quantum tun¬ 
neling through it. 


IV. CONCLUSION 

The calculation of amplification of a weak harmonic 
signal on a noisy background due to the stochastic res¬ 
onance in an RF SQUID loop with a clean ScS con¬ 
tact in the zero-temperature limit shows that the SR in 
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ScS SQUID is possible at any vanishingly small (3l < 1 
because of unusual ’’sharp” potential barrier between 
metastable states whose height is always finite at T = 0. 

With the ScS SQUID temperature rise, the top of the 
barrier smoothes, its height gradually vanishes to zero, 
and the potential of the ScS SQUID becomes much like 
that of’’common” SIS SQUID. Thus, the substantial dif¬ 
ference in classic stochastic dynamics of magnetic flux of 
ScS and SIS SQUIDs emerge at ultra-low temperatures. 
The difference will become more distinctive if quantum 
tunneling contribution is taken into account. These dis¬ 
tinctions should be considered, as useful or ’’parasitic” 
effect, in designing deep-cooled devices involving super¬ 
conducting loops with incorporated Josephson junctions, 
such as SQUIDs and qubits. 
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